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1.  Introduction 


While  the  bulk  of  projectiles  in  use  today  behave  as  rigid  bodies  while  in  flight,  a  notable 
number  of  projectiles  are  purposefully  designed  to  carry  a  liquid  payload.  For  example,  smoke 
screens  delivered  by  artillery  rounds  consist  of  a  typical  spin-stabilized  shell  containing  a 
canister  filled  with  white  phosphorous.  White  phosphorous  is  in  a  liquid  state  when  hot. 

Another  example  is  new  less-than-lethal  projectiles  having  a  concentric  cylindrical  cavity  filled 
with  liquid  that  delivers  this  payload  to  a  target  upon  impact.  Finally,  some  projectiles  are 
designed  to  be  general  payload  delivery  shells,  including  delivery  of  medical  supplies  such  as 
intravenous  fluid  bags. 

Projectiles  with  fluid  payloads  have  been  heavily  researched  over  the  years,  mainly  due  to  the 
fact  that  these  rounds  can  exhibit  severe  flight  instabilities.  Characteristics  of  this  flight 
instability  are  sharp  increases  in  projectile  aerodynamic  angle  of  attack  accompanied  by  large 
changes  in  spin  rate.^ 

Predictions  of  the  instability  induced  by  a  liquid  payload  installed  in  a  projectile  have 
predominantly  been  analyzed  by  linear  fluid  dynamic  theory  subjected  to  a  linear  projectile 
coning  motion.  Unfortunately,  these  well-developed  theories  do  not  directly  mesh  with  standard 
six-degree-of-freedom  (6-DOF)  projectile  flight  dynamic  models.  Hence,  while  a  good  deal  of 
information  is  known  about  the  effect  of  liquid  payloads  on  projectiles,  this  information  has  not 
made  its  way  into  6-DOF  computer  tools.  The  purpose  of  this  report  is  to  bridge  the  gap 
between  the  body  of  literature  on  effects  of  liquid  payloads  on  projectiles  and  6-DOF  projectile 
flight  dynamic  modeling.  The  report  begins  with  a  review  of  rigid  projectile  flight  dynamic 
modeling  along  with  a  description  of  modeling  rotating  liquid  in  a  cavity.  The  two  models  are 
subsequently  integrated  such  that  a  projectile  flight  dynamic  model  with  a  liquid  payload  results. 
This  flight  dynamic  model  is  exercised  on  an  example  shell.  Comparisons  are  made  between  a 
liquid-filled  projectile  and  a  similar  solid  projectile  highlighting  predictive  capability  of  the  new 
model. 


2.  Projectile  Flight  Dynamics 


A  6-DOF  rigid  projectile  model  is  employed  to  predict  the  dynamics  of  a  projectile  in  flight. 
These  equations  assume  a  fiat  Earth.  The  6-DOF  comprises  the  three  translational  components 
describing  the  position  of  the  projectile’s  center  of  mass  and  the  three  Euler  angles  describing  the 
orientation  of  the  projectile  with  respect  to  the  Earth.  Eigures  1  and  2  provide  a  visualization  of 
the  degrees  of  freedom. 


^Karpov,  B.  G.  Experimental  Observations  of  the  Dynamic  Behavior  of  Liquid-Filled  Shell,  BRL-TR-1171;  U.S.  Army 
Ballistics  Research  Laboratory:  Aberdeen  Proving  Ground,  MD,  August  1961. 
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Figure  1.  Projectile  position  coordinate  definitions. 


Figure  2.  Projectile  orientation  definitions. 
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The  equations  of  motion^  for  the  6-DOF  model,  derived  in  the  no-roll  frame,  are  shown  in 
equations  1-4. 
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The  foree  aeting  on  the  projeetile  in  equation  3  comprises  the  weight  force  (W),  aerodynamic 
force,  and  liquid  payload  force  (L).  The  aerodynamic  force  is  split  into  a  standard  (A)  and 
Magnus  (M)  aerodynamic  force.  The  combination  of  forces  is  expressed  in  equation  5. 
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Equation  6  provides  the  expression  for  the  weight  force  in  the  no-roll  coordinate  system. 
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Siegers,  N.;  Kyle,  J.;  Costello,  M.  Nonlinear  Model  Predictive  Control  Technique  for  Unmanned  Air  Vehicles.  Journal  of 
Guidance,  Control,  and  Dynamics  2006,  29  (5),  1179-1188. 
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Equation  7  provides  the  expression  for  the  aerodynamie  foree  in  the  no-roll  eoordinate  system. 
This  foree  aets  upon  the  projeetile  at  the  aerodynamie  eenter  of  pressure. 
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Equation  8  provides  the  expression  for  the  Magnus  foree  in  the  no-roll  eoordinate  system.  The 
Magnus  foree  aets  upon  the  projeetile  at  the  Magnus  force  center  of  pressure. 
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Equations  7  and  8  are  based  on  Mach-number-dependent  coefficients,  the  aerodynamic  angles  of 
attack  given  in  equations  9  and  10,  and  the  total  aerodynamic  velocity  given  in  equation  11. 


d  =  tan-‘(wA/uA). 

(9) 

p  =  tan-‘(vA/u^). 

(10) 

=  VUa  +Va  +Wa  . 

(11) 

The  moment  acting  on  the  projectile  in  equation  4  comprises  the  moment  due  to  the  standard 
aerodynamic  force  (A),  the  moment  due  to  the  Magnus  aerodynamic  force  (M),  the  unsteady 
aerodynamic  moment  (UA),  and  the  liquid  payload  moment  (E)  as  shown  in  equation  12. 
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The  moment  due  to  the  aerodynamic  force  is  expressed  in  equation  13. 
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The  moment  due  to  the  Magnus  force  is  expressed  in  equation  14. 
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The  unsteady  aerodynamic  moments  acting  on  the  projectile  are  expressed  in  equation  15. 
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The  coefficients  used  in  this  model  are  projectile-specific  functions  of  the  Mach  number  of  the 
projectile.  For  fin-stabilized  projectiles,  Magnus  force  and  moment  are  typically  ignored  since 
their  effect  is  rather  small  for  slowing  rolling  projectiles. 

The  dynamic  equations  of  motion  expressed  in  equations  1-15  are  highly  nonlinear.  Due  to  this 
fact,  numerical  integration  is  commonly  used  to  obtain  solutions  to  this  initial  value  problem. 


3.  Liquid  Payload  Moments 


Angular  motion  of  a  projectile  is  altered  by  the  internal  motion  of  the  liquid  payload.  The  liquid 
dynamics  in  this  report  parallel  the  dynamics  given  by  Murphy^  such  that  the  linear  projectile 
theory  gives  a  second-order  differential  equation  for  ^  =  P  -l-  i  a  written  in  terms  of  the  lateral 
force  and  moment 


3 

Murphy,  C.  H.  Angular  Motion  of  a  Spinning  Projectile  With  a  Viscous  Liquid  Payload;  ARBRL-MR-3194;  U.S.  Army 
Ballistics  Research  Laboratory:  Aberdeen  Proving  Ground,  MD,  August  1982.  (See  also  Journal  of  Guidance,  Control,  and 
Dynamics  1983,  6,  280-286.) 
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Solving  this  last  equation  gives  the  linear  sum  of  two  complex  polar  arms: 


=  Ki  e'^‘‘’‘  +  K2 


.iTjPt 


(17) 


where 

ta(K,/K,.)  =  <!,T,P<  . 

Tj-<J/2[1±V1-1/S/, 

Sg  =a"PV4M  if  |H|  «1,|2|«1,  and 
T.H-aS 

g  = _ 1 _ 

'  (2Tj.-a)Tj.P’ 


Consider  a  projectile  containing  a  liquid-fdled  cavity  subjected  to  the  angular  motion  of  a  free- 
flight  projectile.  This  angular  motion  generates  inertial  waves  in  the  contained  liquid,  which,  in 
turn,  impacts  the  projectile  angular  motion  since  liquid  moments  are  now  acting  on  the  projectile 
The  assumption  is  made  that  liquid  moments  are  small  enough  so  the  linear  theory  assumptions 
used  to  obtain  equation  16  are  not  violated. 

X,Y,Z  and  X,  Y,  Z  are  respectively  body  and  nonrolling  coordinate  systems,  where  the  x-axis 
is  the  projectile’s  symmetry  axis.  The  earth-fixed  axes  Xg,  Yg,Zg  has  Xg  initially  along  the 
velocity  vector  and  Z,  Zg  are  initially  downward. 


Moments  due  to  the  liquid  payload  are  assumed  to  be  caused  by  the  coning  motions;  the  linear 
response  of  these  moments  is  modeled  as^4 


Mg  +iNg  =  mg  a 


T,Cg.  K^e’^+T^Cg.  K^e'' 


(18) 


^Cooper,  G.  R.  Moment  Exerted  on  Coning  Projectile  by  a  Spinning  Liquid  in  a  Cylindrical  Cavity  Containing  a  Porous 
Medium-,  ARBRL-MR-3677;  U.S.  Army  Ballistics  Research  Laboratory:  Aberdeen  Proving  Ground,  MD,  June  1988. 
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Clm  will  depend  on  time,  ,  e-,  the  shape  of  the  cavity,  Reynolds  number,  and  the  direction  of 
spin.  Note  that  the  factor(s)  Tj  are  introduced  because  the  liquid  moments  will  vanish  when 
Tj  =  0.  Furthermore,  Clm  causes  rotations  in  the  plane  of  exp^i  (|)j  j  as  well  as  rotations  out  of  the 
plane,  so  the  following  definition  is  introduced; 

^LMj ~  ^LSMj  i  ^LIMj  ’ 

where  Clsm  ^^dCLiM  are  real  and  represent  the  liquid  side  moment  and  the  liquid  in-plane 

-5 

moment  coefficients,  respectively. 

Combining  the  payload  moment  of  equation  19  with  the  aerodynamic  force  and  moment 
produces  a  slightly  more  complicated  differential  equation  for  ^  .  The  projectile  angular 
motion  continues  to  be  an  epicycle  of  equation  16,  which  causes  the  frequency  and  damping  to 
take  the  forms 


T,=T,+T,CL^ 

((f.-l)^CL.„T,f 

fj— l+(mLa  /lx)cLiMj 

£j  =  £aj  +CLSMj  (^L  a"  /Ix )  (2  T  j/a-l)  ' . 


(20) 


The  last  equation  shows  that  Clsm  has  the  same  effect  on  damping  angular  motion  as  does  the 
aerodynamic  damping  moment.  Since  the  coefficient  Clsm  is  usually  positive  for  the  fast  mode, 
this  means  that  Clsm  >  0  and  will  hence  undamp  the  motion.  Similarly,  negative  Clsm  will 
undamp  the  slow  mode. 

In  addition  to  the  two  transverse  liquid  moments  (which  Murphy^  has  shown),  the  axial  liquid 
moment  L  takes  the  form 


L  =  mL  a"  T.  KJ  Cl 


C  =  -C  + 


Ts. 


1- 


3  a" 


(21) 


^Murphy,  C.  H.  Liquid  Payload  Roll  Moment  Induced  by  a  Spinning  and  Coning  Projectile',  ARLBR-TR-0252 1 ;  U.S.  Army 
Ballistics  Research  Laboratory:  Aberdeen  Proving  Ground,  MD,  September  1983. 
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thus  the  roll  moment  is  obtained  direetly  from  the  liquid  side  moment  caleulation. 

When  analyzing  the  impact  of  inertial  waves,  it  is  convenient  to  consider  only  that  part  of  the 
moment  exerted  by  the  liquid  due  to  only  one  of  the  two  projectile  coning  frequencies  Tj .  For 
steady-state  linear  motion,  this  part  will  be  the  total  liquid  moment,  but  for  unsteady  or  nonlinear 
liquid  motion,  this  part  will  be  an  average  contribution  of  the  actual  liquid  moment.  Thus,  the 
total  liquid  moment  will  be  the  sum  of  two  averaged  liquid  moments. 

4.  Liquid  Motion  Model 

The  linear  liquid  analysis  used  for  this  study  is  identical  to  the  analysis  given  by  Murphy  and 
thus  only  a  cursory  review  is  presented  here.  Consider  a  projectile  with  a  cylindrical  payload 
cavity  with  radius,  a  ,  and  height  2c .  The  cylinder’s  axis  is  collinear  with  the  projectile  axis 
with  its  center  located  at  the  projectile’s  center  of  mass.  Linear  theory  is  used  to  predict  the 
liquid  moment  caused  by  coning  motion,  which  has  the  form 

l  =  =KC^  (j  =  l,2),  (22) 

where 

|  =  Kje'^J  =Ke^^ 
s  =  (Sj  +i)Tj  , 

^  =  j)i  ,  and 

K  =  Kj„e'^^»  0  =  1,2). 

The  liquid  motion  is  approximated  quasi-statically  as  the  sum  of  two  steady-state  projectile 
coning  motions.  The  liquid  velocity  components  and  liquid  pressure  have  the  same  dependency 
on  time  and  0  as  the  velocity  components  of  points  on  the  coning  projectile.  Thus,  four  small 
dimensionless  functions,  v ,  w ,  u ,  p  of  r  and  x,  are  introduced; 

V^  =  R(veJ*‘-‘'’)(a^) 

Vg  =  R^we^**"'® j  +  r^/j 
V,=R(ue^"-‘®)(a^^) 

P  =  [R(pej^‘-'®)  +  ^^^rV2](p,ay),  (23) 
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for  the  liquid  cylindrical  velocity  components  V, ,  Vg ,  V^,  and  pressure  P.  These  are  placed  in 
the  linearized  Navier-Stokes  equations  and  the  continuity  equation  to  yield 


(s-i)v-2w  +  a  — =  Re  ' 
^  ^  dr 


a^v  2a^iw 


A.v- 


(s-i)w  +  2v-^  =  Re"' 

(s-i)u  +  a  —  =  Re  '  A^u 
^  ’  dx  " 

5rv  .  5u  ^ 

- iw  +  r —  =  0  , 

dr  dx 


.2  w 

A,w - ^ 


(24) 


where 

d^  d  d^  1  " 

- T  - 1 - T - T  ■ 

dr^  rdr  dx^  r^ 

These  equations  are  solved  by  assuming  the  liquid  has  low  viscosity  so  viscous  effects  can  be 
ignored  except  for  high-gradient  boundary  layers  near  the  cylinder  walls.  This  means  the  inverse 
Reynolds  number  is  taken  to  be  very  small,  Re”'  « 1 ,  so  inertial  wave  solutions  are  modeled  as 
the  sum  of  inviscid,  Re”'  =  0 ,  contributions  plus  viscous  boundary  layer  contributions  of  order 

1/V^. 


v  =  v,+v, 
w  =  w,  +  w^ 

U=Uj+U^ 

P  =  P,+Pv  (25) 


Order  of  magnitude  arguments  applied  to  equation  25  when  substituted  in  equation  24  show  the 
boundary  layer  solutions  near  r  =  a  have  the  following  form: 


w^  =  [^(l  +  is)x/aK- w, 
u^=  (i-s)K  +  Ui 


^(r/a-l)/8. 


^(r/a-l)/8. 


-Jw^dr 


,  5v,  /.  nX 

V, -aoa^^  =  (i-s  —  ,  atr  =  a  . 
5r  ^  ^  a 


(26) 
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Near  the  end  walls  one  will  find 


w.+iu.=-(w,+iu,)e“l'"« 


„P(1+x/c) 


w^-iu^=  Wj -iu; +  2(i  +  s)  — K  e 

d 

Pv  =  0 

^11  r  '' 

Ui  +  5^e-z-^  =  -i(l  +  is)  — Katx  =  ±e  . 

U  X  d 


Note  that  [w^  ,u^  ,p^]  ^  0  exponentially  with  distance  from  the  container  walls;  the  details 
regarding  these  solutions  are  explained  in  an  appendix  in  Murphy.^ 


5.  Liquid  Moments 


The  major  components  of  the  liquid  moment  are  due  to  the  pressure  on  the  walls  of  the  container, 
and  the  viscous  shear  produces  smaller  contributions  to  this  moment.  Therefore,  liquid  moment 
coefficient  is  written  as  the  sum  of  four  terms, 

T  Clm  =  Api  +  fiip,  +  fii,i  +  fii,,  •  (28) 

The  first  two  terms  are  the  pressure  contribution  and  the  remaining  two  terms  are  the  viscous 
shear  contribution.  The  fluctuating  part  of  the  inviscid  P,  which  is  the  primary  contributor  to  the 
inertial  wave  motion,  is  given  by 


Ap=pLaVR  Pi-^K  r{c 


Calculating  the  surface  integral  of  this  with  the  appropriate  moment  arm  over  the  container  wall 
yields  the  moment  contribution  due  to  pressure 


j;j^x[Cpe''<^p+Pv]L^^dx- 


mpi+m  =R  . 

r  dr 

2Ca 
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The  viscous  shear  moment  along  the  lateral  wall  is  calculated  from  the  boundary  layer  solutions 

dx,  (31) 


m,  =  (2KRe)'j:^^ 


■  ai, 

la^^+Cx^— 
dr  dr 


while  the  shear  moment  due  to  the  end  walls  is  given  by 


[2aKRe('j; 

f(w,-iv,)l 

dx 

rdr . 


(32) 


6.  Results 


The  examples  given  here  consider  liquid  payloads  in  the  M483  Army  projectile.  The  liquid 
cavity  is  a  cylinder  with  aspect  ratio  C/a  =  3  completely  filled  with  water.  This  projectile  has  a 
gun  launch  spin  rate,  ^  =  1500/s ,  which  is  large  enough  so  large  Reynolds  number  boundary 
layer  analysis  adequately  governs  the  liquid  physics  along  the  entire  trajectory.  The  range  of  T’s 
for  a  typical  trajectory  of  the  M483  with  a  frozen  (solid)  liquid  payload  is  given  in  figure  3. 


Figure  3.  Coning  frequency  for  the  M483. 

Applying  steady-state  liquid  analysis  to  this  liquid  configuration  reveals  the  liquid  side  moment 
coefficient  for  C/a  =  3,  for  two  Reynolds  numbers  Re  =  5  x  10®,  5x10^  (see  figure  4). 

These  results  depict  typical  Clsm  behavior  where  the  peaks  indicate  large  overturning  liquid 
moments,  possible  resonances,  generated  by  this  contained  liquid  as  a  function  of  the  fast-coning 
Ti  rate. 
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Side  Moment  Resonanee 
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Coiling  Frequency  T 


Figure  4.  Clsm  vs.  coning  frequency. 

This  plot  suggests  that  coning  frequencies  in  the  neighborhood  of  T  =  0.072  will  create 
significantly  larger  side  moments  than  those  found  for  frequencies  outside  this  neighborhood. 
This  result  is  characteristic  of  the  root  cause  of  projectile  instabilities  due  to  liquid  payloads. 

The  aspect  ratio  coupled  with  a  large  enough  Reynolds  number  will  force  inertial  waves  in  the 
coning  fluid  to  generate  large  overturning  liquid  moments.  Increasing  Re  causes  the  liquid  side 
moment  to  increase  rapidly,  resulting  in  a  pronounced  peak  at  a  resonant  frequency. 

The  range  of  frequencies  given  in  figure  4  is  close  to  the  fast  mode  T i  frequencies  frozen  liquid 
found  in  figure  3.  This  is  one  indicator  that  a  flowing  liquid  could  cause  the  fast-coning 
frequencies  to  change  such  that  an  instability  occurs  near  Ti  =  0.072. 

Figures  5  and  6  show  results  of  spin  rate  and  roll  moment  for  the  M483  comparing  a  frozen 

liquid  to  a  flowing  liquid  payload.  The  parameters  selected  for  these  two  plots  are  chosen 

2 

because  they  are  strong  indicators  of  flight  instability  caused  by  liquid  payloads. 

This  example  exhibits  no  flight  instability  caused  by  the  liquid  moments  even  with  the  increased 
magnitude  of  Mx.  The  decreased  roll  rate  and  roll  moment  due  to  the  liquid  is  a  result  of  the 
dominance  of  the  fast  mode  liquid  moment  for  which  Clsm  >  0. 

To  force  flight  an  instability,  take  into  consideration  figure  4  and  let  the  liquid  have  a  decreased 
viscosity  so  that  the  Reynolds  numbers  all  increase  during  flight.  Experience  suggests  that 
undamping  8^  >0.6  means  the  projectile  has  encountered  flight  instability  due  to  liquid 
moments.  An  example  for  the  new  liquid  that  generates  large  overturning  moments  causing  the 
projectile  to  go  unstable  is  now  discussed.  Figure  7  shows  the  time  dependence  of  the  side 
moment  coefficient,  Clsm  i,  rapidly  increasing  neart  1 1  s  . 
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Figure  5.  Comparison  of  projectile  roll  moment  for  frozen  and  flowing  liquid. 


Figure  6.  Comparison  of  projectile  spin  rate  for  frozen  and  flowing  liquid. 
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Figure  7.  Time  history  of  Clsm  i- 

These  results  again  suggest  that  liquid-induced  instability  may  exist  close  to  t  =  1 1  s,  and  the 
cross  plot  shown  in  figure  8  reveals  a  rapid  increase  in  Clsm  i  in  the  neighborhood  of  the  fast 
coning  rate  Ti  =  0.072. 

Similar  results  are  found  for  the  roll  moment  coefficient  Clrm  i,  indicating  that  the  fluid 
resonance  found  in  figures  7  and  8  is  responsible  for  the  large  overturning  moment,  roll  moment, 
and  undamping  exhibited  by  the  fast  frequency  coning  motion.  Figure  9  compares  the  roll  rates 
of  an  unstable  example,  s  >  0.6,  which  occurs  at  time  t  w  1 1  s,  and  a  frozen  liquid  payload. 


Fast  Mode  Clsm 


Fast  Mode  Frequency  Ti 


Figure  8.  Slow  mode  Clsm  fast  coning  frequency  Ti. 
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Figure  9.  Roll  rate  vs.  time  for  flowing  and  frozen  liquid  payloads. 

The  rapid  decrease  in  roll  rate  is  caused  by  the  low-viscosity  liquid,  which  is  forcing  this 
projectile  fast  mode  coning  frequency  to  approach  values  close  to  fast  mode  resonance  near 
Tj  «  0.072.  Figure  10  presents  pitch  rates  of  a  projectile  with  both  frozen  and  flowing  low- 
viscosity  liquid  payloads.  Again,  the  rapid  increase  in  pitch  rate  takes  place  near  time  t  » 1 1  s  in 

the  neighborhood  of  the  fast  mode  resonant  frequency  Ti  =  0.072.  Another  indicator  of  liquid- 
induced  instability  is  the  rapid  increase  in  the  angle  of  attack  given  in  figure  1 1 . 


Figure  10.  Comparison  of  pitch  rates  for  flowing  and  frozen  liquid  payloads. 
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Figure  1 1 .  Angle  of  attack  vs.  time. 


The  data  displayed  in  the  last  five  plots  provide  flight  parameters  that  eorrespond  to  flight 
instability  due  to  the  side  moment  generated  by  the  high  Reynolds  number  liquid  payload.  In 
general,  these  kinds  of  results  occur  when  the  damping  rate  is  8j  >0.6 . 


7.  Conclusion 


Liquid  moments  calculated  from  linear  liquid  motion  undergoing  two-mode  steady-state  coning 
motions  are  used  as  average  moments  when  applied  to  a  nonlinear  6-DOF  time-dependent 
trajectory  model.  These  quasi-static  averages  yield  approximate  predictions  describing  the 
motion  of  projectiles  with  low-viscosity  liquid  payloads  filling  a  cylindrical  cavity.  The 
magnitudes  of  the  liquid  moments  are  dependent  on  the  two  projectile  coning  frequencies,  two 
projectile  coning  undamping  rates,  payload  aspect  ratio,  and  the  liquid  Reynolds  number.  The 
nonlinear  6-DOF  motion  of  a  projectile  is  often  well  approximated  as  the  sum  of  fast  and  slow 
coning  motions.  Thus,  under  these  conditions  the  linear  liquid  analysis  can  predict  when  a 
projectile  exhibits  flight  instability  due  to  a  liquid  payload  by  tracking  the  coning  frequencies 
during  the  6-DOF  numerical  integration  process.  Even  though  such  predications  are  based  on 
steady-state  theory,  they  at  least  serve  as  an  indicator  when  liquid-induced  flight  instabilities  are 
likely  to  exist.  This  hybrid  numerical  tool  provides  reasonable  estimates  of  capturing  the 
physical  effects  of  liquid  moments  on  free-flight  projectiles. 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


a 

c 

f. 

J 

P 

P 

s 


K 


radius  of  the  cylinder  containing  fluid 

half  length  of  the  cylinder  containing  N  subcylinders 

1  +  271  a  Cp/lxCLjjytj 

liquid  mass  hil  =  27ia^  Cp 

projectile  angular  velocity  along  the  projectile  axis 

liquid  pressure  perturbation 
pressure  perturbation 
complex  coning  frequency  s  =  (E  +  i)T 
nondimensionalized  coning  frequency 
gyroscopic  stability  factor 


u 

axial  velocity 

liquid  perturbation  velocity  components 

V 

= 

radial  velocity 

w 

azimuthal  velocity 

Cartesian  unit  vectors 

coning  angle 

gyroscopic  stability  factor 
time 

=  (l  +  i)[2(l  +  w)Re]''^^ 
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Sc 


P 


'cSr, 


2^/l 


+  IS 


1-w 


^  3  + is  ^ 


•  +  I 


+  IS 


IS  J 


growth  rate  per  cycle 

growth  rate  per  cycle  due  to  aerodynamic  moment 
the  dynamic  viscosity  of  liquid 

liquid  mass  density 
variable  in  nonrolling  system 
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